Mon. Not. R. Astron. Soc. 000, 000-000 (2013) Printed 12 April 2013 (MN JAT^K. style file v2.2) 

Mapping UHECRs deflections through the turbulent 
galactic magnetic field with the latest RM data 
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ABSTRACT 

We study the influence of the random part of the Galactic magnetic field on the 
propagation of ultra high-energy cosmic rays. Within very mild approximations about 
the properties of the electron density fluctuations in the Galaxy we are able to derive 
a clear and direct relation between the observed variance of rotation measures and 
the predicted cosmic ray deflections. Remarkably, this is obtained bypassing entirely 
the detailed knowledge of the magnetic properties of the turbulent plasma. Depending 
on the parameters of the electron density spectrum, we can either directly estimate 
the expected deflection, or constrain it from above. Thanks to the latest observational 
data on rotation measures, we build a direction-dependent map of such deflections. 
We find that over most of the sky the random deflections of 40 EeV protons do not 
exceed 1 to 2 degrees, and can be as large as 5 degrees close to the Galactic plane. 
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1 INTRODUCTION 

Magnetic fields play a crucial role in the propagation of 
ultra-high energy cosmic rays (UHECRs) — cosmic rays 
with energies E > 10 19 eV. If the extragalactic magnetic 
fields are not too strong (< 1 nG), the influence of the 
Galactic Magnetic Field (GMF) should be dominant. Var- 
ious observations indicate that the magnetic field of our 
Galaxy c ontain s both regu lar and turbulent components 
(see, e.g., iBeckl (|200ll . 120081 )'). If there were only a regular 
component, a fair knowledge of it would allow one to re- 
store the original positions of the sources from the observed 
arrival directions of CRs. However, the presence of the tur- 
bulent part renders this task very difficult. Unpredictable 
deflections of CRs in the turbulent fields sum up to ran- 
dom displacements that smear the source image. Thus, the 
prospects to reconstruct the actual positions of sources de- 
pend directly on the parameters of the random part of GMF 
(hereafter, rGMF). 

There is quite a substantial number of studies 
concerned with the questi o n of UHECR propagation 
throug h the GMF (IStaneyl dl99 ' 
(|2002f ); IProuza fc Smfdal (|2003h ~ 
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iGiacinti et alj (|2010f ). to name a few). The influence of the 
random component is rather difficult to esti mate and the 
number of papers on this subject is limited (JHarari et al.l 
|2002| ; iTinvakov fc Tkachevl l2005l ; iGolup et al.ll2009l ). 

The rGMF affects not only the deflections of the CRs, 
but also contributes to the scatter of the observed values 
of the Faraday Rotation Measures (RMs) of extra-galactic 
radio-sources. Moreover, the rGMF enters both quantities 
in a very similar fashion: they are proportional to the line- 
of-sight integrals of the perpendicular components of rGMF 
(in the case of deflections) or the parallel component times 
the electron density (in the case of RMs) . Thus, with some 
assumptions about the electron density, one can derive pre- 
dictions/constraints on the random deflections directly from 
the observed RM variances. The aim of this work is to es- 
timate the UHECR deflections due to rGMF using the re- 
cent observational data on the RMs of extragal actic sources, 
namel y the catalogue of RMs obtained by I Taylor et al.1 
(|2009h from the NRAO VLA Sky Survey (NVSS), further 
simply referred to as the NVSS RM catalogue. 

The underlying logic is straightforward. In general, the 
variance of RMs over a sufficiently small patch of the sky 
is composed of the part due to the magnetic fields of our 
Galaxy (this is the part we are interested in) and another 
part due to the intergalactic space, the RM intrinsic to the 
sources, and due to the catalogue errors. The extragalac- 
tic contributions can be isolated and subtracted by mak- 
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ing use of their independence of the Galactic latitude. In 
turn, the Galactic part of the RM variance is a sum of con- 
tributions due to rGMF and the variation of the electron 
density. This relation allows one to either estimate (if the 
second part can be neglected) or put an upper bound on 
the part due to rGMF, which translates directly into an es- 
timate/constraint on the UHECR deflections. The result is 
a position-dependent map of estimates/constraints on the 
UHECR random deflections. Such a map is important in 
the UHECR data analysis, for instance, for a more accu- 
rate prediction of UHECR flux from a given population of 
sources. 

The paper is organised as follows. Section [2] describes 
the data. Section [3] presents our method. The results and 
conclusions are given in section U The closing Appendix 
contains all the technicalities of our method. 



2 DATA 

In this paper we used the NVSS RM catalogue consisting 
of 37,543 RM values of extra-galactic sources. We cleaned 
the catalogue removing the outliers. The algorithm for the 
cleaning procedure was the following: a circle of 3° radius 
was circumscribed around every source and the average RM 
and its variance were calculated for the selected region. If 
the RM of the source was more than two r.m.s. values away 
from the average value, the source was marked as "outlier" . 
In total, 1974 sources were removed after this procedure, 
leaving 35,569 in the reliable set. 

We binned the d ata using the HEALPix pack- 
agqj (|G6rski et alj|2005f ) into a map of resolution A s id e = 
8 in the galactic coordinates with the RING pixel ordering. 
The total number of bins is equal to 768 and the area of 
each bin is 53 sq. deg. This bin size provided the statistics 
of ~ 50 RMs in a bin on average. At the same time, it is 
sufficiently small to neglect the variation of the RM due to 
the change of the coherent magnetic field across the bin. 



3 METHOD 

3.1 Instrinsic variance 

As already mentioned, the observed RM variance a consists 
of two main parts: the variance ctism due to the interstellar 
medium (ISM) we are interested in, and the variance ctqso 
that is caused by other factors, including those intrinsic to 
the sources. We estimated the variance ctqso using the data 
itself by making use of the fact that the latter part should be 
independent of direction, while the latitude (b) dependence 
of ctism should look like 1/sinb due to purely geometrical 
reasons. As we show below, a simple toy model of identical 
cells of turbulent field up to some height z t hr predicts this 
dependence (see section [3~3l . 

In order to estimate ctqso we performed the follow- 
ing calculations. In the Northern hemisphere we fitted 
the observed cr(b) for 36 uniformly spaced values I = 
0°, 10°,. ..,350° with the function /(b) = (A 2 / sin 2 b + 
' 2 treating A and ctqso as free parameters. The bins 



°QSO 



with b < 10° were excluded. We then calculated the mean 
and the r.m.s. values of the resulting 36 values of ctqso- 

The same procedure was independently performed for 
the Southern hemisphere. In both cases we obtained very 
close numbers: ctqso = (12 ± 3) rad m~ 2 . Thus, for the 
following calculations we adopted ctqso = 12 ra d m~ 2 . This 
value is in agreement with other studies, e.g. (|Schnitzelerl 
120101 ) and is largely dominated by the intrinsic errors of the 
NVSS RM catalogue cr N vss = 10.4 ± 0.4 rad m~ 2 . 

Making use of this result, we estimated ctism for each 
bin as follows: 



CTISM 



f 7 ' 
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a > 12 rad m 



(1) 



Both a and ctqso are determined with errors. The error of 
ctqso is cictqso = 3 rad m -2 as explained above. The error of 
a depends on the number n of sources in the bin and equals 
8a = a/y/2n). This implies the error c>ctism in determination 

of CTISM, 



<5ctism = y <5ct qso + 5a 2 



In those bins where the value of ctism inferred from eq. JTJ 
was smaller than <5ctism we have set ctism = <5ctism. This 
includes 29 bins with a < ctqso- 

Finally, we discarded the bins with the number of 
sources less than 10 (130 bins total) located in the "blind 
spot" of the NVSS survey (declination < —40°) and its im- 
mediate vicinity. 

3.2 Rotation measures vs. deflections 

When propagating through a magnetised plasma, the polar- 
isation plane of a linearly polarised electromagnetic wave of 
wavelength A rotates by the angle Aip proportional to the 
square of the wavelength, 



Aip = RM ■ X 2 



(2) 



where the coefficient RM, the rotation measure, can be ex- 
pressed in terms of the properties of the interstellar medium 
(ISM) and permeating magnetic fields. Since the same mag- 
netic field also deflects UHECRs, both processes can be de- 
scribed by similar expressions, functions of the direction r: 



(3) 



RM(f) =ci / dznB 



5i(f) = c 2 



(4) 



where n is the electron density, B\\ and Bi, i — (1,2) 
are the components of the magnetic field parallel and per- 
pendicular to the direction I, respectively, D is the dis- 
tance to the source, and the normalisation constants are 
ci = 0.81 rad cm 3 /(m 2 pc/iG) « 2.7 x 10" 23 rad/^G, c 2 = 
Ze/(E^iG), Z and E being the UHECR charge and energy, 
respectively. For our further estimates we assume that UHE- 
CRs are protons (Z = 1) and have energy equal to 40 EeV, 
thence C2 = 7.5 x 10 -24 rad/(/xG cm). For further reference, 
we also give the ratio 



— = 0.28 cm" 1 . 

Cl 



(5) 



http://hcalpix.jpl.nasa.gov 



These estimates could be easily rescaled for any given com- 
position and energy of UHECRs. 
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The presence of the random component induces vari- 
ances both on the observed RMs and on the deflections of 
UHECRs. Although both mean deflections and mean RMs 
are by definition zero for a random field, their r.m.s. are not. 
The variations in RM are caused by fluctuations both in the 
magnetic field and the electron density, whereas variations 
in S are produced by magnetic fields fluctuations alone. This 
makes the link between the two rather contrived in general. 

Let us separate the regular and random components as 
n — {n)+8n and Bi = {B) i + SBi, and similarly for B3 = B« 
— all these quantities are taken to not depend on the ab- 
solute vertical distance to the observer, see the Appendix. 
Conversely RM and Si already refer to the (square roots of) 
their variances, and there will be no need to prepend a "<5" 
to them. With this notation, and within very mild approx- 
imations we can spell out the connection between (RM 2 ) 
and (S 2 ) = J2{6i6i). 

Schematically, for some positive coefficients a and /3 of 
appropriate units, we show in the Appendix that we can 
write: 



(RM 2 ) = ci (n) 2 (T(8B) 2 ) + a (1(8B) 2 ) + f3 , 



(5 2 ) 



c\ (1(SB) 2 ) 



(6) 



cf. ([3| and (J4}. We also have used the shorthand notation 
I{-k) = f dz-k. In eq. ((6} , the first term comes from the 
variation SB alone, the second depends on both SB and 
Sn (the latter enters through the coefficient a), while the 
third depends only on Sn. The coefficients a and f3 are both, 
therefore, proportional to (Sn 2 ). 

There are several parameters which enter the game at 
this stage, and impact the values of the coefficients a and 
j3. Before dwelling upon these details, we see that there are 
two options. If the properties of Sn are such that both a 
and f3 are small, the corresponding terms can be neglected. 
Then the magnetic fluctuations factorise out altogether and 
we immediately infer 



(S 2 ) = 



d 



c\(n)' 



(RM 2 ). 



(8) 



In this case we link directly (within the limits of the regular 
electron density model) rotation measures with deflections. 
In the opposite case of non-negligible a and /3, there 
are additional, positive by construction, contributions to 
(RM 2 ). Then, the contribution of the magnetic field vari- 
ance (SB 2 ) alone to the RM one would be an overestimation, 
and ([8]) remains valid although only as an upper limit. We 
discuss the goodness of this upper bound below, but first we 
develop a simplified toy model which provides an additional 
intuitive argument to understand the physics behind J8J. 



3.3 A toy random walk model 

The variance of the deflections experienced by UHECRs 
could be estimated within a simple cell model: the random 
field with given strength B r is chopped into identical cells 
with size l ce ii', in each cell the magnetic field is uniform and 
has some random orientation. The electron density n is as- 
sumed to be constant. RMs would then perform a random 
walk along the line of sight, which we can put in formulas 
as 



where fei is a numerical coefficient, N is the total number of 
cells along the line of sight, RMq is the average contribution 
of the elementary cell — one can think of this as © with 
fdz — > /ecu- 

The deflections in this model are described by a very 
similar formula: 



(10) 



5 = k 2 S VN 



where k 2 is another coefficient, and So is the deflection in 
the elementary cell (once again, take ((4} and replace J dz — ¥ 

^coll). 

The exact value of the coefficients k\ , k 2 were calculated 
numerically. First we simulated the random walk for the RM 
accumulation. Taking into account that only the parallel 
component of the magnetic field contributes proportionally 
to cos (f>i, where 4>i is the angle between the random direction 
of the magnetic field and the direction of the line of sight, 
then, fci reads 



fei 



2V 

Y, cos< 

1=1 



N 



(11) 



where (f>i is a random variable in [0, n] range. The number 
of steps N was taken equal to 100,000 and we averaged over 
10,000 realisations, obtaining finally 

kt = 0.56. 

The coefficient k 2 could be obtained in the very same 
fashion, but now there is a two-dimensional random walk 
process and also it is the orthogonal component of the mag- 
netic field that should be taken into account. One finds 



\( N 
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/ N 


2 


4/1 J2 sin fa cos & 


) + 


Y sin 4>i sin ipi 


) 


V V=i 




\! = 1 


/ 



k 2 = J-i^-: — - — — , (12) 



N 

where <f>i is again a random variable in the [0, n] range, 
whereas ipi has domain [0, 27r]. We obtain 

k 2 = 0.63. 

Putting everything together, the relation between the 
RM and the deflection along the same path is 



5 ft; 0.31 cm 



,1 k 2 RM 
ki (n) 



0.35 cm" 1 ^-. (13) 

(n) 



RM = kiRMoVN, 



(9) 



This result is to be compared with the more accurate expres- 
sion J8} where, for E — 40 EeV and Z = 1 (protons), the 
coefficient turns out to be cijd w 0.28 cm -1 , cf. eq. ([5}. The 
two expressions are therefore in close agreement, the more 
so in light of the observational and theoretical errors we are 
faced with: the toy model correctly distillates the physics, 
and provides clear and intuitive insight for our idea. 



3.4 General case 

We turn now to discussing the parameters which enter in 
a and /3, and the interpretation of |[SJ|. For concreteness we 
focus on the simplest, and mos t corroborated by the data, 
scenario (|Armstrong et al.lll995l ; lHaverkorn et al.ll20081 ). We 
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thus employ a Kolmogorov electron density spectrum to ob- 
tain 

{RM 2 } ~ c\ {n) 2 {X(SB) 2 ) (14) 

+ c\ [3.1C 2 l 2 /3 (l(SB) 2 ) + l.W 2 n Dl 5 /3 {B} 2 ] , 

where C 2 is the amplitude of the Kolmogorov spectrum, 
and the momentum IR cutoff is given by qo — 2tt/Iq, with lo 
the largest scale of turbulence (with a given spectral index); 
finally, D is the path length along which the electron den- 
sity turbulence extends (and, as we expect for the magnetic 
field one, it will roughly follow a 1/ sin(fe) distribution). The 
details of the calculation are given in the Appendix. 

As previously mentioned, there are two possibilities. If 
the magnetic field variance dominates this expression, then 
there is a simple direct relation {RM 2 } oc (T(8B) 2 ). The 
second option is that the {B} 2 term wins, in which case the 
measured RMs would be due to a combination of fluctua- 
tions in the electron density and the magnetic field, which 
we can not disentangle; the contribution of the magnetic 
field variance {SB 2 } to the UHECR deflections would then 
be smaller than in eq. ([8|, so the latter would be an over- 
estimation and would, therefore, provide one with an upper 
bound on the deflections. 

We can also turn this argument around: assuming a 
given set of parameters which describe the electron density, 
its fluctuations and the average magnetic field, we can ask 
which would be the contribution of the last term (indepen- 
dent of the unknown {SB 2 }) to the measured RM. If this 
is actually larger than what is observed, the conclusion is 
that the set of parameters used is not consistent with the 
RM data. Still, eq. © gives a valid upper bound for a given 
electron density. 

Consider the parameters that enter eq. (|14[l . First, we 
fix the value of the path length D that should be in- 
serted into our calculations. The variances get contributions 
only in regions where considerable magnetic fields and elec- 
tron density are present. Both quantities are usually as- 
sumed to decay very rapidly wi th the height h above the 
Galactic pl ane: like sech 2 (fe/d) (jCordes fc Lazid 12002! ) ; or 
exp(— h/d) (jGaensler et alj2008r ). albeit with possibly some- 
what different vertical scales d. To describe the electron den- 
sity d istribution we use the NE2001 model (jCordes fc Lazid 
[2002J) modified with an increased vertical scale (n(h — 0) w 
0.014 cm" 3 , d = 1.8 kpc), which is implied by recen t ob- 
servations JGaensler et alj 120081 : IPshirkov et al.ll201ll l. We 
adopt this vertical scale and cut the integrals at the value 
D = d/sin(6). Notice that since all three terms in eq. (|14p 
are roughly proportional to D, their relative magnitude does 
not depend strongly on the adopted value. 

Moreover, notice that we do not treat explicitly the h- 
dependence of the average (n) itself, which would also be- 
have as sech 2 (h/d) or exp(— h/d), although we do keep the 
dependence with direction. In doing so we are making an 
approximation equivalent to substituting the steep suppres- 
sion with a sharper step function localised at h = d, with 
the advantage that we can calculate everything analytically. 
This boosts the deflection we obtain by at most around 30%, 
which is not large an error in view of the accuracies we have 
at hand. 

The mean magnetic field (B) entering eq. (|14|l can 
be inferred from the existing models of the coher- 



ent G alactic field. We ad opt (B) estimated from the 
model (|Pshirkov et al.l 2011T); alternativ ely, one could use, 
e.g., the model ( Jansson fc Farrarll2012l ). 



The statistical properties of Sn are no t very 
well known. The an al yses of lArmstrong et al.l (|l995T ) 
and Irlaverkorn et al.1 (|2008l ) have obtained C 2 = 
10~ 3 m~ 20/3 and an IR scale of l = 100 pc (at large b). 
Again, these values are expected to change with direction 
because of the morphology of the galaxy, but the data are 
uncertain and we will rely on the simplifying approximation 
that these parameters remain the same across the sky. With 
these values, for a typical average magnetic field of 1 [iG 
and b = 90°, the third term in eq. (|14[) amounts to about 17 
rad/m 2 , which overshoots the observed RM in a large part 
of the sky, although not by much. We have therefore a sit- 
uation described above when the parameters are not quite 
consistent with the observed RMs. If this inconsistency is at- 
tributed to the numerical uncertainties, one concludes that 
the third term in eq. (|14[) dominates and thus eq. (|8]) is an 
overestimate. 

Let us now look at slightly different parameters of the 
electron density fluctuations. So far we have always assumed 
exact Kolmogorov turbulence all the way up to the IR cut- 
off lo of 100 pc. Nevertheless, the observational scenario is 
far from being clear, and there is the possibility that the 
spectrum actually flattens at larger scales; there are i ndica- 
tions of this flattening for the turbulent magnetic field (jRegisJ 
[2011J) , although it is uncertain what the link between the two 
is at large scales. Within our formalism it is immediate to 
estimate its effects: as a concrete case, taking the spectral 
index a = 3 from 1 pc to 100 pc suppresses the contribution 
of the third term in eq. (|14jl to the total RMs by a factor 
of 100" 1 , bringing it down to about 4 rad/m 2 , which is 
below the intrinsic threshold we adopt. In this case |(HJ be- 
comes a much better estimation of the deflections: if we ask 
that the third term in eq. (|14|l contributes at most to 20% 
of the observed RM, then any RM above about 20 rad/m 2 
would provide a direct connection to {S 2 }. 

The last comment is on the second term in eq. (|14[) . 
This term has the same dependence on {SB 2 } as the first 
one, so the two can in principle be combined together, which 
would change (make smaller) the coefficient in eq. ([8]). Al- 
though this would strengthen the constraints on the deflec- 
tions for a given value of RM, it would also make them 
dependent on the (poorly known) properties of the electron 
density fluctuations. For this reason we do not include this 
term i n our final constraints. Usin g again the fiduc i al val - 
ues of lArmstrong et al.1 <|l995h and Irlaverkorn et al.1 (|2008T ) 
we find 1.7QC n l ~ 0.080cm~ 3 , which is larger than the 
model average of NE2001; the latter is, (n) « 0.012cm" 3 ). 
Using instead the spectral index a — 3 from 1 pc to 100 pc 
suppresses this contribution down to 0.017 cm -3 , thereby 
bringing it around the mean value of NE2001. 

To summarise, within present uncertainties of the pa- 
rameters of the electron density fluctuations, both third and 
second terms in eq. (|14|l can be sizeable, but their values 
cannot be reliably determined at present. Because of these 
uncertainties, we base our constraints on eq. (|8]) which we 
treat as an upper bound on deflections {S 2 }. 
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Figure 1. The map of displacements ^/(9 2 ) of protons with en- 
ergy E = 4x 10 19 eV due to the random component of the GMF 
as follows from cq. JSJ. This map should be considered as an up- 
per bound and may turn to an estimate if the contribution of the 
electron density fluctuations are small, as explained in the text. 



4 RESULTS AND CONCLUSIONS 

As explained above, eq. (JSJ) provides one with an upper 
bound (or an estimate, if the contributions of the fluc- 
tuations of the electron density can be neglected) on the 
UHECR deflections (S 2 ) in the random component of the 
Galactic magnetic field. These deflections can be translated 
into the experimentally observable displacements {8 2 } of the 
UHECR arrival directions at the Earth with respect to the 
direction to the true position of the source. Th e relation 
between the two quantities is (jHarari et alj|2002h : 

(e 2 } = (5 2 )/3. (15) 

In Fig. [T] we show the map of deflections \J {8 2 ) of pro- 
tons with energy E — Ax 10 19 eV due to the random com- 
ponent of the GMF calculated according to eq. @ where we 
have used the value C2/C1 = 0.28 cm -1 . This map should be 
considered an upper bound, which may turn into an estimate 
if the contribution of the electron density fluctuations will 
shown to be small in the future. The displacements \J{6 2 ) 
are bound to 1 to 2 degrees for the majority of directions, 
although closer to the galactic plane the bounds are relaxed 
to ~ 5°. This map is the main result of our paper. 

One can see from Fig. [l] that the limit on deflections 
does not depend on the Galactic longitude in any regu- 
lar way. In particular, the direction towards the Galactic 
center does not look any special. On the contrary, there is 
a regular dependence on the Galactic latitude b. This de- 
pendence is fitted well with the phenomenological function 
^/(sin 2 (&) + B), with the constants A = 0.5, B = 0.13, 
although the fluctuations are large. For the purpose of de- 
scribing the upper bound, the constraints of Fig. [T] can be 
majorated by the analytic function 



l c 



sin 2 (6)+0.15 



(16) 



everywhere except several bins. 

The resul ts obtained in the paper are i n good agreement 
with results of lTinvakov fc Tkachevl (|2005r ) . In the latter pa- 
per the deflections in the turbulent field were derived using 
properties of the rGMF inferred from few observations com- 
ing from different sets of data. It was shown that for 40 
EeV protons such deflections are small, not exceeding 1.5° 
for the better part of the sky. The most important differ- 



ence with that work, where one was obliged to use a specific 
model (i.e., spectrum) for the rGMF fluctuations to draw 
some quantitative conclusion, is that here we were able to 
place robust upper limits from a single, thereby consistent, 
set of data. 

In conclusion, a word of warning is in order. While the 
constraints shown in Fig. [l] are conservative in the sense 
that we have treated all the uncertainties in a way that in- 
creases the deflections (and thus weakens the constraints), 
they only apply to the contribution of the region of ±2 kpc 
around the Galactic disk where the electron density is large 
and constrained by observations. One cannot exclude, by 
the arguments presented in the paper, that there might exist 
more extended regions around the Galaxy that carry small 
electron density and substantial magnetic fields such that 
their contributions to the RMs are small, but the contribu- 
tion to deflecti ons of UHECRs i s larger than in eq. ((Sj like 
in the model of lAhn et al.l (|1999f ). Constraints on such mod- 
els require a different approach, which we leave for future 
study. 
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APPENDIX 

We describe here in details the formal steps that read to the 
relation (JSJ between the variances of RM and deflections. 



Definitions 

The two quantities we want to calculate are 



RM(f) = a / dznB 



Si(f) — C2 dze,j Bj ; 

Jo 

more specifically, we need the variances of these quantities. 
In principle, both the random variations in electron density 
and magnetic field are present. Thus, we split n = (n) + Sn 
and Bi = {B) i + SBi, where i — (1,2) denotes the two 
components perpendicular to the line of sight, and similarly 
for S3 = By • 

In what follows we approximate the actual distance 
(from the observer) dependence of such quantities is that of 
a pure step function. This is certainly a good approximation 
as it captures its key feature: a sharp drop beyond a given 
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distance D. Indeed, we do expect that both the electron den- 
sity and magnetic field (regular components and fluctuations 
alike) follow the smoothed distribution of i onised gas in the 
galax y, which is ty pically coded as sech 2 (:r) IjCordes fc Laziol 
120021 ) or exp(-a;) (|Gaensler et alj2008h (a; is the vertical dis- 
tance from the observer normalised to d); both these func- 
tions do not vary exageratedly up to x ~ 1 and rapidly decay 
beyond that point — the mistake we are making in treating 
these as step functions is of, at most, around 30%, which is 
well below our other uncertainties. We are therefore allowed 
to average, along the line of sight, the regular components 
(n) and (B) within the integrals and bring them outside. 
This also implies that the spectra of fluctuations will not 
depend on the absolute position, but only on the scale. 
Thus, we define 



Zy = {8B i (ro)6B j (r + r)) 



<Tqe 



ifq PB (g) 

2qi 



QiQj 



where we used the Fourier transform 



SBi(f) 



i3 irq 

a qe 



SBi(q), 



and the spectrum is defined as 
(5BM6B*(P)) = ^0- (Sa q -f ) 8*(q p) . 
Similar definitions apply to Sn: 
8n(f) —Id qe rT ^8n{q) , 



and 



Pn(q) s3/ = 



(5n(q)5n* (p)> = — —8 {q-p). 

Notice that ly oc oy since otherwise it is zero (parity) . Writ- 
ing the integrals explicitly, they all reduce to the prototype 

I / J J QJ 2 ~irq cos u> P B (OJ I -, ?i 

h = / d cos (pdfidqq e w -y II 

where we call z the line of sight axis. The angular integrals 
can be solved to give 

f do _ ,. /sin or cosgr sin or \ 

h = h = 2tt / -ift? H - + -=± 3^- , 

J q \ qr q^r z q°r° J 

T T , f dq „ , . fcosqr s'mqr 
h = J|| =4tt / —P B {q) ' 



Our assumptions on the properties of these fluctuations with 
vertical distance amount to assuming that the spectra we 
have hitherto defined do not depend on such quantity, but 
only on the (modulus of the) correlation length r. 

The next step is to integrate along the line of sight: 

f D 
(8i8i) — c| / dzdz'eije ik (8Bj(z)8B k {z')) 

Jo 

rD ru 

= 2c^ / du / dr(/i + h) 

Jo Jo 

= AtycI [ —P B (q) (qDSi{qD) + cos qD 



sin qD 



here we have defined r = z' — z and u = z' + z, and Si(x) is 
the sine-integral function. 



A simplified result 

A very similar result can be obtained for the RM provided 
we ignore the fluctuations in the electron number density. 
Let us first write the full expression for the variance of RM, 
containing all pieces: 



(RM 2 ) =c\ dzdz 



{n) 2 (8B\\{z)8B\\{z')) 



(17) 



+ (8n{z)8n{z')){B)\ + ((SB ]l Sn)(z)(8B ]] 8n)(z')) 



where the first term, {SB 2 }, is the contribution of rGMF 
proper, the second one, {Sn 2 }, is that of the electron density 
fluctuations alone, and the last one represents the convolved 
variations. 

If the variance on the electron density fluctuations is 
small, then the only significant term left in this lengthy ex- 
pression is the first, leading to 



(RM 2 )=c 2 {n) 2 / dzdz'(5B {l (z)5B ll (z')) 
Jo 

= 47r Cl j ^P fl (g) (qDSi(qD) + cos qD - 2 + ^|p 



To further simplify this expression, one can expand for large 
distances and small coherence length qD 3> 1 to obtain 

(S 2 )=4nc 2 2 J^P B (q)^ q D + . 
{RM 2 } = ^c\(n) 2 J ^P B (q) (|«U 

This demonstrates now the two quantities are, in first ap- 
proximation, directly related, and one can infer the deflec- 
tions caused to the UHECRs without knowing anything about 
the spectrum of the magnetic field. Of course this result 
makes use of the assumption that there are not any large 
electron density fluctuations, the goodness of which we dis- 
cuss in the next section. 



The convolution 

We now want to determine the contribution of the convolu- 
tion 



{(SB ll Sn)(z)(6B ll 5n){z')} 



-- I d 3 kd 3 k'e- 1 ^~ l¥i ' 



{{5B\\5n)(k)(8B\\5n){k')) 



= {A*y 



qe 

dp dq 



i(q+ P )(z'-z) P«{p) PBJq) I ^ _2l 

" g 2 



P n (p)P B (q) 
p q pr 



p° 2q 3 
sinpr j svaqr 



cos qr 



This has to be integrated twice over z and z , ox r and u as 
usual. One reasonable assumption, which simplifies greatly 
the (euphemistically) unwieldy result, is that the integral is 
going to be dominated by the IR cutoffs, since observations 
tell us that the spectra are almost certainly red. The integral 
can then be split in two, for p <g q and p ^> q depending 
on which IR cutoff is the most important one. The result 
can (now) be expanded for large qD if this is the dominant 
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cutoff, or, conversely for large pD to find 



(RM 2 ) C (2tt) 3 Ci 



dp dq 
P Q 



Pn(p)P B (q) 



IS. 

3 q 

D_ 
P 



p< q, 
p> q- 



Compare these with the integrated {8n 2 } and (SB 2 ) in 
the same limit of large pD and qD — recall that I(*) = 
Jdz*: 



(l(Sn) 2 ) 



-&-K 



I 



p p 

D 



{I(6BY) = -2-k* J : fP B (q) q 

one can write the convolution part as either something pro- 
portional to one or the other. In any case, the two limits 
differ only by a factor of 2/3, so we can take either one: the 
mistake we make is at most this much. Physically, the first 
option luight be more justified, that is, turbulence in elec- 
tron density extends beyond, at larger scales, that for the 
magnetic field, so, for concreteness, we choose this in what 
follows. The main result then reads 

(l(8nSB) 2 ) = -^L(T{5Bf) f ^P n (p) . 

As this relation illustrates, everything can be predicted with 
the mild assumption that we know the features of the elec- 
tron density spectrum, and there is no necessity to specify 
the properties of the magnetic field random fluctuations. No- 
tice that thus far the only approximation we have employed 
is that the IR dominates the fluctuation spectra. 



Results 

We need to know the average electron density (n) and reg- 
ular magnetic field (B), and the spectrum of the electron 
density fluctuations P n (q) , in order to be able to extract the 
deflections from the RM data bypassing the magnetic field 
fluctuations. 

The spectrum P n can be written as 
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{Sn 2 )(f) 



C^nC 



{q 2 +q 2 Y> 2 



where C^ is the amplitude of a Kolmogorov spectrum for the 
electron density fluctuations, and the momentum IR cutoff is 
given by qo = 2n/lo, with lo the largest scale of turbulence 
(with a given spectral index). This integral can be solved 
exactly, and it gives 



{Sn 2 )(r) 



2^ 2 C 2 n ( r \ 
r(a/2) \2qo) 



(a-3)/2 



K, 



(a-3)/2 



(gor) 



These definitions enable us to compute all the terms in (1171) . 
As an exa mple, following t he re sults of lArmstrong et al.l 
(|l995T ) and lHaverkorn et al.l (120081 ) . for Kolmogorov electron 
density fluctuations (a = 11/3), we find 



(RM 2 ) 



{ [(nf + 3.1C^ 2/3 ] (1(SB) 2 } + l.lC 2 n Dll /3 (B} 2 } 



which is the result (|14[) reported in the main text, and jus- 
tifies the use of ((6j). The ensuing discussion then follows. 



